Helix or Coil? Fate of a Melting Heteropolymer 
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We determine the probability that a partially melted heteropolymer at the melting temperature 
will either melt completely or return to a helix state. This system is equivalent to the splitting 
probability for a diffusing particle on a finite interval that moves according to the Sinai model. 
When the initial fraction of melted polymer is /, the melting probability fluctuates between different 
realizations of monomer sequences on the polymer. For a fixed value of /, the melting probability 
distribution changes from unimodal to a bimodal as the strength of the disorder is increased. 

PACS numbers: 82.35.-x, 02.50.Cw, 05.40.-a 



Thermally-induced helix-coil transitions in biological 
heteropolymers exhibit intriguing thermodynamic and 
kinetic features [l|, [2|, |3, |4[ that stem from, for example, 
different intrinsic melting temperatures of the monomeric 
constituents, say, A and B, and quenched random distri- 
butions of monomers along the chain. Due to this ran- 
domness, an arbitrarily long A-rich helix region with a 
high local melting temperature can act as a barrier to 
hinder the melting of the entire chain into a random coil. 
As the temperature is raised, the nucleation of multiple 
coils may occur at distinct points along the condensed 
helical chain. Nevertheless, as considered in Ref. |5j, the 
simplified situation in which there is a single coil portion 
starting at a free end of the chain (Fig. Q]) highlights the 
essential physical role of randomness on melting kinetics. 
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FIG. 1: Illustration of a partially melted heteropolymer with 
the interface between the melt (left) and helix (right) at Y. 

Our primary result is that the probability that a par- 
tially melted heteropolymer returns to its native heli- 
cal state (and the complementary probability that the 
chain denatures) exhibits huge sample-specific fluctua- 
tions. Consequently, the average and the typical behav- 
iors are completely different, and neither is representative 
of the behavior of a single chain. Our analysis is based on 
making an analogy [5| between the motion of the bound- 
ary point (Fig.[T]) in a partially melted heteropolymer and 
a particle diffusing in the presence of random force — the 
Sinai model [6(. 

Let E+ (Y ) be the melting probability for a heteropoly- 
mer of (contour) length L with a partially melted segment 
of length Y, while E- (Y) = 1 - E+(Y) is the probability 
that this heteropolymer condenses. The melting proba- 



bility is equivalent to the splitting probability E + {Y) that 
the boundary point starts at Y and reaches L without 
ever reaching (and correspondingly for E-(Y)). For 
a random potential that corresponds to a heteropoly- 
mer, the splitting probabilities are different for each se- 
quence of monomers. Moreover, the resulting distribu- 
tion of splitting probabilities changes from a single peak 
at its most probable value for weak disorder (with a 
delta- function peak in the absence of disorder), to dou- 
ble peaked with most probable values close to and 1 for 
strong disorder. Thus much care is needed to interpret 
experimental data of heteropolymer melting kinetics. 

Although the monomer sequence in a heteropolymer 
is typically correlated over a finite range, we make the 
simplifying assumption that such correlations are absent. 
The location of the equivalent boundary point between 
the coil and helix portions of the chain then diffuses in 
the interval [0, L] in presence of a quenched and random 
position-dependent force F(x), with 



(F(x)) = F 0) 
(F(x)F(x'))-F$ = o* F 6(x- 



x'). 



(1) 



The mean force Fq > when the temperature T exceeds 
the heteropolymer melting temperature T m (and Fq < 
for T < T m ), while o~f is proportional to the difference 
in melting temperatures for A and B homopolymers [5(. 
Following the analogy with splitting probabilities, the 
probability that a homopolymer (no random potential) 
ultimately melts when a length Y is initially in a random 
coil state is pj 
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where, for later convenience, we define z = e~ ' D and 
/ = Y/L. Here D is the diffusion coefficient of the bound- 
ary point between the melt and helix (Fig. []}, and the 
velocity v > for T > T m and v < for T < T m . For 
T = T m , corresponding to v = 0, these splitting proba- 
bilities simplify to E+{Y) = 1 — E-(Y) = /, the classical 



result for unbiased diffusion [7j. We now focus on the 
corresponding behavior for a heteropolymer at the melt- 
ing temperature where the boundary point moves in the 
random potential defined by Eq. (fTJ). 

We first present an optimal fluctuation method [8j 
to describe the statistical features of V(E + (Y)) to de- 
scribe the ultimate fate of a single chain. The basis of 
this approximation is to replace each realization of the 
Sinai model by an effective environment whose mean bias 
matches that of the Sinai model (Fig. [21). Then the dis- 
tribution of bias velocities in the continuum limit is 



7 is the viscosity coefficient. We then use the fluctuation- 
dissipation relation D = kT/j to rewrite Eq. ([6]) as 



V(E + ) = - 
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■ exp 
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with a = a F L/[8(kT) 2 } a dimensionless measure of the 
strength of the Sinai potential relative to thermal fluctua- 
tions. The important feature of this splitting probability 
distribution is its change from unimodal to bimodal as 
the disorder parameter a increases past a c = 1 (Fig. [3]) . 
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(3) 



Since there is a one-to-one connection between the bias in 
the effective model and the splitting probability E+(Y), 
we may now convert the distribution of bias velocities to 
the distribution of splitting probabilities by 



V(E+) = P(v) 



dv 
dE~ 



(4) 




FIG. 2: A typical realization of the potential in the continuum 
limit of the Sinai model (smooth curve) and the potential in 
the corresponding effective medium (dashed). 

To perform this transformation, we use Eq. @ to solve 
for v as a function of E + for fixed /. This inversion is 
feasible only for / = | , 1 , 1 , j\ , because solving a poly- 
nomial up to quartic order is involved. For the simplest 
case of a particle starting at the midpoint, / = |, then 
@ becomes E+ = 1/(1 + z). Inverting gives 
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Substituting these results into Eq. Q and using (j4|) gives, 



V{E+) = 



2D 



y^na^LE+E. 
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We now relate the velocity variance a 2 to the disorder 
in the Sinai model. For Sinai disorder, the mean-square 
potential difference between two points separated by a 
distance L is (Fig. \2§ 



U, 



eff 



(F(x')F(x"))dx'dx" 



2 r 



Thus there is a net force F c s ~ U c g/L ~ ap/yL, from 
which we infer the velocity scale a v ~ <?f j (■yvL) , where 
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FIG. 3: (color online) Optimal fluctuation prediction for the 
splitting probability distribution for the Sinai chain when the 
particle starts at Y — L/2 (top) and Y — L/3 (bottom, loga- 
rithmic scale). 

In the critical case a = a c the maximum at E+ = 1/2 
is quartic and the distribution is close to uniform for 
0.2 < E+ < 0.8. Thus any value of E + in this range 
is nearly equally probable. For a > a c , two maxima 
emerge continuously at E + = \ [1 ± \Ji{a — a c )] . Curi- 
ously, although the disorder-average value (E+) = 1/2, 
for a > a c the probability distribution has a minimum 
at E + = 1/2. In this strong disorder limit, the underly- 
ing mechanism for the bimodality of the distribution is 



that a typical realization of the environment has a net 
bias (see Fig. [2]) so that the splitting probability is cither 
close to zero or to one. 

For / = i, a similar calculation gives the splitting 
probability distribution shown in Fig. O For weak disor- 
der, the splitting probability is peaked at a point close to 
the value ^ that arises in the case of no disorder. When 
a is sufficiently large, the splitting probability has peaks 
near and 1, corresponding to strong disorder where a 
typical configuration has an overall net bias. The inter- 
mediate regime of a gives the strange situation in which 
there is a peak at one extremum but not the other. 

For a general starting point, the inversion v = v(E + ) 
and thus the form of V(E + ) can be determined asymp- 
totically in the limits \vL/D\ — > oo. For vL/D — > — oo, 
Z — ► +oo, and Eq. @ becomes E + ~ z^( 1 ^/)// '. Solving 
for v(E + ) and -j^- and using these results in Eq. (|4]), the 
splitting probability distribution starting from Y = f L 
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As a is varied, this distribution changes from having a 
peak at the most probable value of E + to a peak at 
E+ = 0. In the complementary limit vL/D — ► +oo, 
corresponding to z —* and E+ — > 1, the distribution is 
again {JSJ> , but with / — ► 1 — / and £ + — ► 1 — J5+. 

We now give an exact solution for the splitting prob- 
abilities that qualitatively substantiates the heuristic re- 
sults given above. Consider the Fokker-Planck for the 
probability distribution for a particle that moves in the 
Sinai potential, P t = D[P XX — /3(F(x)P) x ], where j3 = 
1/kT and P(x,t)dx is the probability that the particle 
is in the range [x, x + dx] at time t and F(x) is the force 
at x. Suppose that particles are injected at a constant 
rate at Y to maintain a fixed concentration P(Y) = 1 
when absorbing boundary conditions at x = and x = L 
are imposed. Then the stationary solution to the Fokker- 
Planck equation is 



P(x) = < 



J Y exp[— /3J-y(x)] dx 
f Y exp[— [5!Fy{x)] dx 

J^exp[-/3J r y(x)]da; 
f Y exp[— /3JV(x)] dx 



x<Y, 

(9) 
where Ty{x) = J Y dx'F(x'). 

From the steady current J = —D{P X — /3(F(x)P)], the 
splitting probabilities E±(Y) are simply 



E. 



(10) 



where t± = l/J±, with J± being the flux to the bound- 
aries at and L, respectively. Thus r± define the "resis- 
tance" of a finite interval with respect to passage across 



the two boundaries and are given explicitly by 



-#Fo(x) dx ^ 



T+ 



L-Y 



e~ 0:Fo{x) dx. (11) 



We emphasize that the force integrals !Fq are independent 
variables in each integrand so that r_ and r + are also 
independent random variables. 

The quantities defined in (fTT|) are the continuous-space 
counterparts of Kesten variables [9(] that play an im- 
portant role in many stochastic processes. For exam- 
ple, negative moments of r describe the positive mo- 
ments of a steady diffusive current in a finite Sinai chain 
|10L I 111 . Il2l Il3 . Il4l | . A surprising feature is that the 
disorder-average currents J± scale with system length L 
as L -1 / 2 and are much larger than Fickian which scale 
as L _1 in homogeneous environments. Thus Sinai chains 
have an anomalously high conductance despite the fact 
that diffusion is logarithmically confined [6( and positive 
moments of r grow as e ^ [HJ, HE H ES Ei | ■ 

Let ^-(t-) and ^ + (t + ) denote the distribution func- 
tions of the random variables r_ and r+, respectively. 
Then the moment generating function of the splitting 
probability E + can be written as: 
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(12) 

Integrating over dr+ , we formally change the integration 
variable from r + to E + to give 
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where r = E-/E+. From this expression we may read 
off the following form of the probability density V(E + ) 
of the splitting probability E + : 



V(E + ) = -L r 'dr_T_*_(r_)*+(rr_). 

E+ Jo 



(13) 



For a random potential without any bias, the distribu- 
tion ^-(t-) is given by (see, e.g., Ref. [lJl): 
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where a_ = /3 2 cr|,/L/4 and chw, denotes the hyperbolic 
cosine. Then ^ + (r+) is obtained from (fT4]) by replacing 
r_ — > t + and a_ — > a + = /3 2 cr|,(l - f)L/A. 

Substituting *_(r_) and *+(r + ) into Eq. (JT3J) and 



integrating over dr_ , we find 
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After cumbersome but straightforward manipulations, 
one of the integrals can be performed to recast the split- 
ting probability distribution as: 
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where rj = sh _ [(ch u)/^/r\. For the symmetric initial 
starting point, a + = a_ = a, the distribution changes 
from unimodal to bimodal as a increases beyond a c w 
0.815 (Fig. UK . 




FIG. 4: (color online) Exact distribution V{E + ), Eq. (T16J1 . 
for a- — q+ = 1/4 (solid curve) a± — 0.815 (dashed), and 
a± = 1.2 (dash-dotted). Thin curves are the corresponding 
asymptotic results of Eq. (|17[) . 

To extract the asymptotic behavior of T J (E+) in the 
limits E+ — > 1 and E+ — > note that because the in- 
tegrand in Eq. (fl"6|) contains an oscillating cosine term, 
only the behavior near u = should matter. Thus as- 
suming chu w 1 and expanding ?y 2 in a Taylor series in 
u up to second order, we find the following asymptotic 
representation for E+ — ► 1: 
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where z + = (1 + \/ ' E + ) / \J ' E-. This asymptotic form 
agrees quite well with the exact result in Eq. (jTSJ) , not 



only when E + — > 1, but also for moderate values of E + 
(Fig. [4]). Similarly, we obtain the asymptotics of V(E + ) 
for E + — * merely by interchanging all subscripts ± to 
=F in (|17[) . and z_ is obtained from z + by interchanging 
the subscripts ± to =F in the latter. When E + ^ 1 or 
1 — E + <C 1 respectively, Eq. (fTT)) reduces to 
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Our optimal fluctuation result ([7|) resembles the above 
form apart from logarithmic and numerical factors. 

In summary, by making an analogy to first passage on 
a Sinai chain, we find that the evolution of a partially 
melted random heteropolymer at the melting tempera- 
ture is controlled by the sequencing of monomers. Each 
heteropolymer realization has a unique kinetics and final 
fate that is not representative of the average behavior 
of an ensemble of such polymers. A related lack of self 
averaging was recently found in anomalous diffusion [l5j . 
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